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Abstract: We perform a global analysis of cosmological observables in generalized cos- 
mologies which depart from ACDM models by allowing non-vanishing curvature fi^ ^ 0, 
dark energy with equation of state with w / -1, the presence of additional relativistic 
degrees of freedom AiV re i, and neutrino masses £l v ^ 0. By combining the data from 
cosmic microwave background (CMB) experiments (in particular the latest results from 
WMAP-7), the present day Hubble constant (HO) measurement, the high-redshift Type-I 
supernovae (SN) results and the information from large scale structure (LSS) surveys, we 
determine the parameters in the 10-dimensional parameter space for such models. We 
present the results from the analysis when the full shape information from the LSS matter 
power spectrum (LSSPS) is included versus when only the corresponding distance mea- 
surement from the baryon acoustic oscillations (BAO) is accounted for. We compare the 
bounds on the neutrino mass scale in these generalized scenarios with those obtained for 
the 6+1 parameter analysis in ACDM + m v models and we also study the dependence of 
those on the set of observables included in the analysis. Finally we combine these results 
with the information on neutrino mass differences and mixing from the global analysis of 
neutrino oscillation experiments and derive the presently allowed ranges for the two labo- 
ratory probes of the absolute scale of neutrino mass: the effective electron neutrino mass 
in single beta decay and the effective Majorana neutrino mass in neutrinoless /3j3 decay. 
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1. Introduction 

It is now an established fact that neutrinos are massive and leptonic flavors are not sym- 
metries of Nature [1,2]. In the last decade this picture has become fully proved thanks 
to the upcoming of a set of precise experiments. In particular, the results obtained with 
solar [3-12] and atmospheric neutrinos [13, 14] have been confirmed in experiments using 
terrestrial beams: neutrinos produced in nuclear reactors [15, 16] and accelerators [17-20] 
facilities have been detected at distances of the order of hundreds of kilometers [21]. 

The minimum joint description of all the neutrino data requires mixing among all the 
three known neutrinos (f e , i^, z/ T ), which can be expressed as quantum superpositions of 
three massive states Vi {i = 1, 2, 3) with masses m^. This implies the presence of a leptonic 
mixing matrix in the weak charged current interactions [22,23] which can be parametrized 
as: 



(I \ / ci3 s 13 e- iS c P \ I C2i si2 0\ (e^ 0\ 
U = c 23 S23 ■ 1 • -ai 2 ci 2 • e** , (1.1) 

\0 -S23 c 23 / \si3e i5cp ci3 / \ 1/ \ 1/ 



where Cij = cos 6^- and Sij = sin Oij. In addition to the Dirac-type phase Sep, analogous 
to that of the quark sector, there are two physical phases rji associated to the Majorana 
character of neutrinos and which are not relevant for neutrino oscillations [24,25]. 

Given the observed hierarchy between the solar and atmospheric mass-squared split- 
tings there are two possible non-equivalent orderings for the mass eigenvalues, which are 
conventionally chosen as 



As it is customary we refer to the first option, Eq. (1.2), as the normal (N) scheme, and 
to the second one, Eq. (1.3) , as the inverted (I) scheme; in this form they correspond 
to the two possible choices of the sign of Am^. In this convention the angles 0^ can be 



vn\ < 7712 < ^3 with Am^j (Am| 2 — Am^) with (Am^ > 0) ; 
TO3 < mi < m2 with Am^ <C |Am 31 ~ Am| 2 | with (Am 31 < 0) . 



(1.2) 
(1.3) 
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taken without loss of generality to lie in the first quadrant, 6ij £ [0, 7r/2], and the phases 
5 C p, Vi G [0,2tt]. 

Within this context, Am^, [AmJ-J, 9\2, and 6*23 are relatively well determined from 
oscillation experiments [26-29], while only an upper bound is derived for the mixing angle 
#13 and barely nothing is known on the phases and on the sign of Am^. Furthermore 
neutrino oscillation data provides as unique information on the absolute neutrino mass 
scale a lower bound 



Y j m>\J\£*r% 1 \ forN (1.4) 

^m J >2^|Am| 1 | fori (1.5) 

(1.6) 

Conversely the neutrino mass scale is constrained in laboratory experiments searching 
for its kinematic effects in Tritium f3 decay which are sensitive to the so-called effective 
electron neutrino mass [30-32] 

2 _ 2|T7|2 22 2,22 2,2 2 /-, m 

m v e = 2^, m i\ U ei\ = c i3 c i2™i + c 13 s 12 m 2 + s 13 m 3 , (1.7) 

i 

At present the most precise determination from the Mainz [33] and Troitsk [34] experiments 
give no indication in favor of m Ue / and one sets an upper limit 

m Ue < 2.2 eV , (1.8) 

at 95% confidence level (CL). A new experimental project, KATRIN [35], is under con- 
struction with an estimated sensitivity limit: m Ue ~ 0.2 eV. 

Direct information on neutrino masses can also be obtained from neutrinoless double 
beta decay (0v/3/3) searches provided they are Majorana particles. In the absence of other 
sources of lepton number violation in the low energy lagrangian, the 0i/(3(3 decay amplitude 
is proportional to the effective Majorana mass of u e , m ee , 



m PP = 



c\ z c\ 2 m x + c\ z s\ 2mi + s 2 13 e~ iS ^ 



(1.9) 



which, in addition to the masses and mixing parameters that affect the tritium beta decay 
spectrum, depends also on the phases in the leptonic mixing matrix. The strongest bound 
from 0vf3f3 decay was imposed by the Heidelberg-Moscow group [36] 

m ee < 0.26 (0.34) eV at 68% (90%) CL, (1.10) 

which holds for a given prediction of the nuclear matrix element. However, there are large 
uncertainties in those predictions which may considerably weaken the bound [37] . A series 
of new experiments is planned with sensitivity of up to m ee ~ 0.01 eV [38]. 

Neutrinos, like any other particles, contribute to the total energy density of the Uni- 
verse. Furthermore within what we presently know of their masses, the three Standard 
Model (SM) neutrinos are relativistic through most of the evolution of the Universe and 
they are very weakly interacting which means that they decoupled early in cosmic history. 
Depending on their exact masses they can impact the CMB spectra, in particular by al- 
tering the value of the redshift for matter-radiation equality. More importantly, their free 
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streaming suppresses the growth of structures on scales smaller than the horizon at the 
time when they become non-relativistic and therefore affects the matter power spectrum 
which is probed from surveys of the LSS distribution (see [39] for a detailed review of 
cosmological effects of neutrino mass). 

Within their present precision, cosmological observations are sensitive to neutrinos 
only via their contribution to the energy density in our Universe, VL v h 2 (where h is the 
Hubble constant normalized to Hq = 100 km s" 1 Mpc -1 ). Q, u h 2 is related to the total 
mass in the form of neutrinos 



Therefore cosmological data mostly gives information on the sum of the neutrino masses 
and has very little to say on their mixing structure and on the ordering of the mass states 
(see Ref. [40] for a recent update on the sensitivity of future cosmological observations to 
the mass ordering.) 

There is a growing literature on the information extracted from cosmological obser- 
vations on the neutrino mass scale, starting with the analysis performed by the different 
experimental collaborations [41-45]. The basic observation is that, besides variations due 
to the observables considered, the bounds on the neutrino mass obtained depend on the 
assumptions made on the history of the cosmic expansion, or in other words, on how many 
parameters besides the ACDM model are allowed to vary when analyzing the cosmological 
data. Additionally depending on those assumptions some observables need or not to be 
considered in order to account for degeneracies among the parameters (for some recent 
analysis see [46-48]). 

In this article we present the results of a global analysis of cosmological observables in 
owCDM + AAr re i + m u cosmologies which depart from ACDM models by allowing, besides 
neutrino masses Q. v ^ 0, non-vanishing curvature 7^ 0, dark energy with equation of 
state with together with the presence of new particle physics whose effect on the 

present cosmological observations can be parametrized in terms of additional relativistic 
degrees of freedom AiV re i. In particular this extends the most general analysis of Ref. [48] 
by accounting also for non-flatness effects. We adopt a purely phenomenological approach 
in analyzing the effect of a non-vanishing spatial curvature without addressing its origin. 
However it is worth mentioning that, within inflationary models which produce the simple 
initial conditions here considered it is difficult to end up with a significant f2/% [49]. We 
describe in Sec. 2 the different cosmological observables included in these 10-parameter 
analysis as well as our statistical treatment of those. The results of the analysis are pre- 
sented in Sec. 3 where we discuss the differences obtained when the full shape information 
from the LSS matter power spectrum is included versus when only the corresponding dis- 
tance measurement from the baryon acoustic oscillations is accounted for. We also compare 
the bounds on the neutrino mass scale in these owCDM + AA^ re i + m„ scenarios with those 
obtained for the 6+1 parameter analysis in ACDM + m u models and we also study the 
dependence of those on the set of observables included in the analysis. These results are 
combined with the information on neutrino mass differences and mixing from the global 
analysis of neutrino oscillation experiments in Sec. 4 to derive the presently allowed ranges 
for the two laboratory probes of the absolute scale of neutrino mass: the effective neutrino 
mass in single beta decay m Ve and the effective Majorana neutrino mass in neutrinoless /3/3 
decay m ee . We summarize our conclusions in Sec. 5. 




(1.11) 
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2. Cosmological Inputs and Data Analysis 



Parameter 


symbol 


Hubble Constant 


Ho 


Baryon density 


n b h 2 


Dark matter density 


Q c h 2 


Scalar spectral index 


n s 


Optical Depth at Reonization 


T 


Amplitude of scalar power spectrum at k = 0.05 Mpc -1 


A s 


Total neutrino mass 




Dark energy equation of state parameter 


i=l,3 
u 


Effective number of extra relativistic degrees of freedom 


AN Ic i 


Spatial curvature density 





Table 1: Cosmological parameters used in our most general analysis, h = H /100. We denote the 
cosmology characterized by these parameters as owCDM + AiV ro i + m v . 

We consider cosmologies owCDM + AiV re i + m v characterized by the free parameters 
listed in Table 1. All parameters are as usually defined in the literature with the exception 
of Ai^i. Our definition of extra relativistic degrees of freedom accounts for the fact that 
we have evidence of the existence of three and only three standard neutrino species which 
mix due to mass oscillations [21]. Their contribution to the energy budget of the universe 
is included in £l u h 2 . A7V re i parametrizes the contribution of additional relativistic massless 
states of any spin to the radiation energy density. For convenience that contribution is nor- 
malized to the one from a spin 1/2 weakly interacting massless state. That normalization 

is AiV rc i. 

In these cosmologies several parameter degeneracies appear in any of the cosmological 
observables. First, any experiment that measures the angular diameter or luminosity dis- 
tance to a single redshift is not able to constrain because the distance depends not only 
on Ofc, but also on the expansion history of the universe. Thus for a universe containing 
matter and vacuum energy, one needs to combine at least two absolute distance indicators, 
or the expansion rates out to different redshifts to break this degeneracy. Furthermore 
when dark energy is dynamical, co ^ — 1 , a third distance indicator is required. Finally 
the presence of extra relativistic degrees of freedom AN TC \ changes the matter-radiation 
equality epoch, a change that can be compensated by the corresponding modification of 
the matter density Q m h 2 . As a result, AiV re i and Q m h 2 are strongly degenerate unless a 
fourth distance indicator provides us with an independent constraint on £l m h?. 

In our analysis we include the results from the 7-year data of WMAP [41] on the tem- 
perature and polarization anisotropics in the form of the temperature (TT), E-mode po- 
larization (EE), B-mode polarization (BB), and temperature-polarization cross-correlation 
(TE) power spectra for which we use the likelihood function as provided by the collabora- 
tion 1 . A number of CMB experiments have probed smaller angular scales than WMAP. 
In particular we consider the results from the temperature power spectra from the Cosmic 

x We notice that, although the models considered do not generate any B-mode polarization, in order to 
account for the information from EE and low-1 TE data, the BB power spectrum must also be included in 
the analysis because WMAP provides the combined likelihood for the low-1 TE, EE and BB spectra [50] 
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Background Imager (CBI) [51], the Very Small Array (VSA) [52], BOOMERANG [53] and 
the Arcminute Cosmology Bolometer Array Receiver (ACBAR), [54]. In order to avoid 
redundancies among the CMB data sets, we follow the procedure in Refs. [42,55]. We use 
seven band powers for CBI (in the range 948 < I < 1739), five for VSA (894 < t < 1407), 
seven for BOOMERANG (924 < t < 1370), and sixteen band powers of ACBAR in the 
range 900 < t < 2000. We do not include the results of the Background Imaging of Cosmic 
Extragalactic Polarization (BICEP) [56] experiment whose bands overlap excessively with 
WMAP and from QUaD [57] which observes the same region of sky as ACBAR and it is 
less precise [55]. Furthermore we do not include in the analysis the polarization results 
of these experiments. As mentioned above, in the analysis of WMAP we use the likeli- 
hood function as provided by the collaboration. For the other CMB experiments we build 
the corresponding likelihood functions from the data, covariance matrix and window func- 
tions given by each experiment. We compute theoretical CMB predictions using the fast 
Boltzmann code CAMB [58,59]. Following the procedure outlined in Ref. [43] whenever 
it is required we account for the Sunyaev-Zel'dovich (SZ) effect by marginalizing over the 
amplitude of the SZ contribution parametrized by the model of Ref. [60]. We assume a 
uniform prior on the amplitude as < Asz < 2. 

We also include the results from Ref. [61] on the present-day Hubble constant, Ho = 
74.2 ± 3.6 km s _1 Mpc -1 where the quoted error includes both statistical and systematic 
errors. This measurement of Hq is obtained from the magnitude-redshift relation of 240 
low-z Type la supernovae at z < 0.1. We include this result as a Gaussian prior and neglect 
the slight cosmology dependence [47] of this constraint. 

The results from luminosity measurements of high-z Type la supernovae are included 
as presented in the compilation of the supernova data called the "Constitution" sample in 
Ref. [62] which consists of 397 supernovae and it is an extension of the "Union" sample 
[63]. With these data we build the corresponding likelihood function without including 
systematic errors whose precise values are still under debate [62]. In our analysis we 
marginalize over the absolute magnitude of the supernovae with a uniform prior. 

Finally we also include the results from the matter power spectrum as derived from 
large scale structure surveys in two different forms. In one case we use the measurement 
of BAO scale obtained from the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and 
the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) [44]. In the other we include 
the full power spectrum of the SDSS DR7 survey [45] (which we label LSSPS). 

For the analysis including the BAO scale, we use as input data the two distance ratios 
at z = 0.2 and z = 0.35 presented in Ref. [44] and build the corresponding likelihood 
function using the covariance matrix as given in that reference. As discussed in Ref. [44] 
the distance ratios can be considered as measurements of d z = rs{zd) / 'Dy(z) and apply to 
any of the considered models. r s (zd) is the comoving sound horizon at the baryon drag 
epoch and Dy{z) = [(1 + z) 2 D 2 A cz/H (z)] 1//3 with Da is the angular diameter distance 
and H(z) is the Hubble parameter. However in their fitting procedure the value of d z 
is obtained by first assuming some fiducial cosmology, extracting the value of Dy(z) and 
then computing rs(zd)/Dy(z) with r s (zd) evaluated by that fiducial cosmology using the 
approximated formula of Eisenstein & Hu [64] for Zd- As discussed in Ref. [48] this approx- 
imated formula is not strictly valid for the extended cosmologies which we are considering. 
We correct for this effect by a) exactly evaluating the redshift at baryon drag epoch by 
using Eq.(B.5) in Ref. [48] in the extended cosmologies, and b) correcting for the use of 
the approximate formula in the presentation of the data by rescaling the predictions by a 
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factor rj d a pp rox )/rf d (zd exact) (we prefer to rescale the predictions since the covariance 
matrix is given for the data as presented). 

In our second analysis we include the full SDSS DR7 data which consists of 45 bins, 
covering wavenumbers from h m 'm ~ 0.02 /^Mpc ^ to & max = 0.2 /iMpc -1 (where k m i n and 
&max denote the wavenumber at which the window functions of the first and last data point 
have their maximum). In this analysis we use the likelihood function as provided by the 
experiment. Together with the linear matter power spectrum it requires a smooth version 
of it with the baryon oscillations removed. We construct such no-wiggle spectrum for the 
extended cosmologies here considered from the linear matter power spectrum computed 
by CAMB using the method based on the discrete spectral analysis of the power spectrum 
described in Appendix A.l of Ref. [48]. 

We also perform comparative analysis including only the ACDM parameters plus mas- 
sive neutrinos (first seven in Table 1) fixing uj = — 1 AN Te \ = £7^ = for different combi- 
nations of the above observables. 

Additional constraints on the cosmological parameters can be obtained if one includes 
in the analysis information on the growth of structure from other low redshift data. Among 
others the small scale primordial spectrum determined from Lyman-alpha forest clouds or 
the priors on the amplitude of mass fluctuations derived from different galaxy cluster sam- 
ples. We have conservatively chosen not to include those in our analysis because generically 
these results are subject to model dependence assumptions which render them not directly 
applicable for the most general cosmologies here consider. 

With the data from the different samples included in a given analysis and the theoreti- 
cal predictions for them in terms of the relevant parameters x, we construct the correspond- 
ing combined likelihood function. In Bayesian statistics our knowledge of x is summarized 
by the posterior probability distribution function (p.d.f.) 

P[Xl >- JC(D\x')7r(x'\V)dx'- 1 > 

it{x\V) is the prior probability density for the parameters. In our analysis we assume a 
uniform prior probability for the x parameters in Table 1. For ^2 m v and AiV re i we imposed 
that they should be both > 0. Following standard techniques in order to reconstruct 
the posterior p.d.f. Eq. (2.1) we have developed a Markov Chain Monte Carlo (MCMC) 
generator which employs the Metropolis-Hasting algorithm including the adapting for the 
kernel function to increase the efficiency Full details are given in Appendix B of Ref. [65] . 
For each combination of data we generate O(50) chains in parallel and verify its convergence 
by studying the variation of the Gelman- Rubin i?-parameter [66] imposing as convergence 
criteria R - 1 < 5 x 10~ 2 . 

3. Results of the Cosmological Fits 

Our results for the two analysis in owCDM + AiV re i + m u cosmologies are presented in 
Figs. 1- 2 and in Table 2. In Fig. 1 we show the marginalized one-dimensional probability 
distributions for the ten independent parameters obtained from Eq.(2.1) as 2 



Pi-dim(^) = J dx k ^ip(x\D,V) . (3.1) 



2 Technically this is obtained from the MCMC chain by discretizing the parameter space and counting 
the fraction of points in each cell. 
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owCDM+AN rel +m„ 
CMB+HO+SN+LSSPS 

owCDM+AN re ,+ m„ 
CMB+HO+SN+BAO 



2.5 5 2 

AN rel Em, (eV) 



Figure 1: Constraints from our global analysis on the cosmological parameters of ou;CDM+AiV re i+ 
m u for the analysis including CMB+HO+SN+LSSPS (solid red) and CMB+HO+SN+BAO (dotted 
blue). The different panels show the marginalized one-dimensional probability distributions for all 
parameters. For the neutrino mass see also Fig. 3. 



For convenience we show the information on the normalization of the scalar power spectrum 
in terms of the derived as parameter which parametrizes the expected root mean square 
amplitude of the matter fluctuations in spheres of radius R = 8/i _1 Mpc. 

The best fit values given in the second and fourth columns in Table 2 are those for 
which Pi_dim( 2; i )est ) is maximum. The allowed ranges at a given CL, xf^ in < x\ < xp^ ax , 
are obtained from the condition 



% , max 



CL [£j jm i n < Xi < Xj >max ] — / Pl-dim(xi) (3.2) 



CL 



*^i,min 



with Pl-dim(^n) = Pl-dim(^maJ ( 3 - 3 ) 



or ■ 

KJ1 z,mm 



for Xi = AN rel ,^2m u (3.4) 



where (3.4) is used when there is no solution for condition (3.3). 

Equivalently we define the marginalized two-dimensional probability distribution func- 
tions 

P2-dim(Xi,Xj) = J dxk^ijpixlD,^) , (3.5) 

and from these, we obtain the two-dimensional credibility regions with a given CL as the 
region with smallest area and with CL integral posterior probability. In practice they 
are obtained as the regions surrounded by a two-dimensional isoprobability contour which 
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CMB+HO+SN+BAO 


CMB+HO+SN+LSS-PS 




best 


la 


95% CL 


best 


lcr 


95% CL 


H km/s/Mpc 

n b h 2 x 100 
n c h 2 

' v o 
T 


76.2 
2.205 
0.131 
0.961 
0.086 
0.787 


+3.0 
-2.8 

+0.057 
-0.050 

+0.018 
-0.013 

+0.021 
-0.015 

+0.011 
-0.015 

+0.091 
-0.073 


+5.7 
-5.6 

+0.103 
-0.105 

+0.036 
-0.023 

+0.040 
-0.030 

+0.026 
-0.028 

+0.135 
-0.179 


74.4 
2.239 
0.128 
0.971 
0.083 
0.824 


+2.8 
-2.9 

+0.059 
-0.046 

+0.024 
-0.009 

+0.019 
-0.017 

+0.016 
-0.011 

+0.051 
-0.048 


+5.6 
-5.6 

+0.095 
-0.108 

+0.042 
-0.018 

+0.037 
-0.033 

+0.030 
-0.023 

+0.097 
-0.105 




-0.006 


+0.010 
-0.009 


-0.022 < n k < 0.016 


-0.011 


+0.008 
-0.009 


-0.028 <n k < 0.007 


UJ 


-1.17 


+0.19 
-0.21 


-0.62 < uj + 1 < 0.18 


-1.12 


+0.21 
-0.20 


-0.57 < uj + 1 < 0.26 


AiV rc l 


1.2 


+1.1 
-0.61 


0.08 < AiV rc i < 3.2 


1.3 


+1.4 
-0.54 


0.21 < AiV rc i < 3.6 


J2 m u (eV) 




< 0.77 


< 1.5 




< 0.37 


< 0.76 



Table 2: Constraints from our global analysis for owCDM + A iV ro i +to„ cosmologies. We show the 
values for the best fit parameters and the corresponding lcr (68%) and 2cr (95%) allowed intervals. 



contains the point of highest posterior probability and within which the integral posterior 
probability is CL. We plot in Fig. 2 the 68% and 95% CL two-dimensional credibility regions 
for the last four parameters in Table 1. For the analysis including CMB+HO+SN+BAO 
(full regions) and including CMB+H0+SN+LSSPS (void regions). 

Because of the degeneracies present in these cosmologies one finds, as expected, a 
degradation in the constraints of the standard parameters (ie those of the ACDM model) 
when compared with the analysis performed within the ACDM priors for the same set of 
observables (see for example table 1 in Ref. [41]). As seen in the figure this particularly 
affects the determination of Q c (or equivalently Q m ) as a consequence of the well-known 
degeneracy in the predictions of the CMB spectra between AJV rc i and Q m . This is so 
because a simultaneous change of both can leave untouched the redshift for matter-radiation 
equality, which is well constrained by the ratio of height of the third and first peaks in the 
CMB spectra. This degeneracy is broken by the addition of the HO prior as well as the 
independent determination of £l m from the distance information from LSS, either using 
only BAO or the full power spectrum. It is interesting to notice that we find that data 
is better described by allowing a non-zero amount of extra radiation even though it is 
only at most a 2a effect. This implies, for example, that models with extra light sterile 
neutrinos are favoured by the data (as discussed in Ref. [67] in the context of flat cosmologies 
with a cosmological constant) even for these owCDM models. Most conservatively we can 
read the results as a 2a upper bound on AA rc i < 3.2 (3.6) for the analysis including 
CMB+HO+SN+BAO (CMB+H0+SN+LSSPS). 

We also find a widening in the allowed range of ns as a consequence of its degeneracy 
with AArei and with the dark energy equation of state uj. Both change the ratio of generated 
power at low versus large angular scales in the CMB spectrum, an effect that can be offset 
by a change in the spectral index. Conversely we find that in these cosmologies uj is 
considerably less constrained than in ljCDM scenarios for which a 95% range —0.089 < 
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In\ (eV) AN rel co Q k 




Void regions CMB+HO+SN + LSSPS 



Figure 2: Constraints from our global analysis for owCDM + AiV re i + m v cosmologies for 
the analysis including CMB+HO+SN+BAO (full regions) and for for the analysis including 
CMB+HO+SN+LSSPS (void regions). We show the 68% and 95% CL two-dimensional credibility 
regions for the last four parameters in Table 1 

uj + 1 < 0.12 is obtained from the analysis of CMB+BAO+SN [41]. 

The normalization of the power spectrum as parametrized by as is mostly affected 
by the presence of neutrino masses. Their main effect is to reduce the amplitude on the 
power spectrum on free-streaming scales therefore decreasing erg. There is also a residual 
degeneracy between eg and Qk mostly associated with the fact that allowing for non- 
flat cosmologies permits to increase the amount of dark matter in the form of neutrinos 
without affecting £l c and therefore minimizing their indirect impact on the CMB spectra. 
As a consequence we find that there is a correlation between the allowed range of neutrino 
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mass and f2fc as it is seen in the corresponding panel of Fig. 2. This leads to a somewhat 
wider allowed range of ^ when compared to the results obtained from the analysis of 
CMB+BAO+SN in owCDM scenarios -0.019 < n k < 0.0072 [41]. 

Figures 1 and Fig. 2 also display clearly the differences in the results obtained when 
the full shape information from the LSS matter power spectrum is included versus when 
only the corresponding distance measurement from BAO is accounted for. We see that, 
with the expected exception of the neutrino mass (and correspondingly of as), both sets 
of data lead to comparable precision on the determination of the cosmological parameters. 
Concerning the neutrino masses we find that neither of the two analysis show any evidence 
for neutrino mass and the best fit point is obtained for ^ m u = 0. However the 95 % 
upper bound obtained when using BAO, ^2,ra v < 1.5, is tighten by a about a factor 2, 
J2 m v ^ 0.76, by considering instead the full LSSPS. 




1.5 2 

I m, (eV) 



Figure 3: Constraint on £m„ as a function of the CL for the different analysis as labeled in the 
figure. 

We plot in Fig. 3 the bound on ^ m u for the two analysis in the owCDM + AiV re i + m v 
cosmologies at a given CL together with the corresponding results from different anal- 
ysis performed in the framework of ACDM + m u models. The corresponding 95% CL 
bounds are listed in Table. 3. We find that for the same combination of observables 
CMB+HO+SN+BAO (CMB+HO+SN+LSSPS) the bound for a ACDM + m v scenario is 
m u < 0.61 (]P m u < 0.35) which is a factor ~ 3 (2) tighter than the corresponding one 
obtained in owCDM + AiV re i + m u cosmologies. However, we also find that at lower CL 
ACDM + m v scenarios are better fitted with a non vanishing m u when the information 
from CMB (and HO) is combined with the information from LSS surveys. This is, however, 
at most a la effect associated with the slight mismatches between the best fit values of 
the cosmological parameters obtained in the analysis of the observables in the context of 
the ACDM. This illustrates the well-known fact that overconstrained scenarios are more 
"sensitive" to small fluctuations in the data, due either to statistics or to an optimistic esti- 
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Model 


Observables 


Era, (eV) 95% Bound 


owCDM + AiV rc i + m v 


CMB+HO+SN+BAO 


< 1.5 


ouCDM + AiV re i + m v 


CMB+HO+SN+LSSPS 


< 0.76 


ACDM + m„ 


CMB+HO+SN+BAO 


< 0.61 


ACDM + m„ 


CMB+HO+SN+LSSPS 


< 0.36 


ACDM + m„ 


CMB (+SN) 


< 1.2 


ACDM + m„ 


CMB+BAO 


< 0.75 


ACDM + m„ 


CMB+LSSPS 


< 0.55 


ACDM + m„ 


CMB+HO 


< 0.45 



Table 3: 95 % upper bound on the sum of the neutrino masses from the different cosmological 
analysis. The analysis within ACDM + m v models including only CMB data or in combination 
with SN yield the same 95% bound. 

mate of the systematic uncertainties. Consequently, even if more conservative, the bounds 
derived on more general scenarios are more robust against these effects. 

4. Combination with Oscillation Data 

We present in this section the allowed ranges for the sum of the neutrino masses and the 
two laboratory probes of the absolute scale of neutrino mass: the effective neutrino mass 
in single beta decay m Ve and the effective Majorana neutrino mass in neutrinoless f3/3 
decay m ee , obtained from the combination of the cosmological analysis discussed above, 
with the information from the global analysis of solar, atmospheric, reactor and accelerator 
longbaseline (LBL) neutrino experiments in terms of flavour oscillations between the three 
neutrinos [26]. 

Our starting point is the \ 2 function from the oscillation analysis 

Xo( A "4l> Am 31> 012, 6*13, #23, ^CP) = XSolar+KamLAND( A ™21, #12, #13) + AcHOOz( A ™31> #13) 

+XATM+LBL( Am L Am L #12, #13, #23, <$Cp) I 4 - 1 ) 

=>• Xo{mve,m ee ,^m n ) (4.2) 

where the last step is obtained after marginalization over Ara^j and #23 and allowing for 
variation of the two phases r/i and 772 within their full range. 

In Fig. 4 we plot the 95% allowed regions (for 2 dof) in the planes {m Ve ^m v ) and 
(m ee ,Y2rn u ) as obtained from the marginalization of Xo( m u e , m-ee, with respect to 

the undisplayed parameter in each plot. In the figure we also show superimposed the single 
parameter 95% bounds on rn v i from the different cosmological analysis described in the 
previous section. The figure illustrates the well-known fact that currently for either mass 
ordering the results from neutrino oscillation experiments imply a lower bound on m Ue . On 
the contrary m ee is only bounded from below for the case of the normal ordering while full 
cancellation due to the unknown Majorana phases is still allowed for the inverted ordering. 

In order to obtain the global combined ranges we first define a one parameter equivalent 
Xc(J2 m i/) function [68] for a given cosmological analysis from the condition that it leads 
to the same CL intervals than the corresponding marginalized one-dimensional probability 
distribution function: 







Cosmo+ Oscillations 








95% Ranges 




Model 


Observables 


m„ e (eV) 


m ee (eV) 


Sm„ (eV) 


owCDM 
+AiV rcl + m v 


PMR-I-HO-I-SN-I-R AO 


N [0.0047- 0.51] 
1 [11.1)4/ — U.OlJ 


N [0.00-0.51] 
1 [U.U14 — U.OlJ 


N [0.056- 1.5] 
1 [U.Uyo — 1.0] 


oujCBM 
+AN Td + m„ 




N [0.0047 - 0.27] 
I [0.047 - 0.27] 


N [0.00 - 0.25] 
I [0.014 - 0.25] 


N [0.056 - 0.75] 
I [0.098 - 0.76] 


ACDM 4- m 


CIMR-I-Hn+SN+RAO 


N [0.0047 - 0.20] 
I [0.048 - 0.21] 


N [0.00 - 0.20] 
I [0.014 - 0.21] 


N [0.056 - 0.61] 
I [0.097- 0.61] 


ACDM -I- m 


("'MR-J-Hfl-l-^N-l-T 


N [0.0047- 0.12] 


N [0.00-0.12] 


N [0.056 - 0.36] 


I [0.047-0.12] 


I [0.014 - 0.12] 


I [0.098 - 0.36] 


ACDM + m, 


CMB f+SN") 


N [0.0047 - 0.40] 
I [0.047 - 0.40] 


N [0.00 - 0.40] 
I [0.014 - 0.41] 


N [0.056- 1.2] 
I [0.098 - 1.2] 


ACDM + m v 


CMR+RAO 


IN [U.UUOZ — U.ZOJ 

I [0.047 - 0.25] 


IN [U.UU — U.ZOJ 

I [0.014 - 0.25] 


IN [U.UOD — U. i OJ 

I [0.099 - 0.75] 


ACDM + ra. 


CMB+LSSPS 


N [0.0047- 0.18] 
I [0.048 -0.19] 


N [0.00-0.18] 
I [0.014 - 0.19] 


N [0.056 - 0.55] 
I [0.099 - 0.55] 


ACDM + m v 


CMB+HO 


N [0.0047-0.14] 
I [0.047-0.16] 


N [0.00-0.14] 
I [0.014 - 0.16] 


N [0.056 - 0.44] 
I [0.097 - 0.45] 



Table 4: 95% allowed ranges for the different probes of the absolute neutrino mass scale from the 
global analysis of the cosmological data with with the results from oscillation experiments. The 
analysis within ACDM + m v models including only CMB data or in combination with SN yield the 
same 95% ranges. 



where CL is obtained from Eq. (3.2) with Y m Vi = ^f^ax anc ^ x ?min = u when the lower 
bound for that CL is (which imples that the function Xc(Y m v) ^ s sm gl e valuated). 
When xf^ in / the function Xc(Y m v) takes the same value for J2 m ^i = x ?max an d 



Finally we construct 



Y m Vi — x i,mw. 



Xo+c("V> TO ee,X^J = Xo{rn Ue ,m ee ^m Ui ) + XcC^rn u ) , (4.4) 

from which we obtain the the 2a 1-dim allowed ranges for m Ve , m ee , and Y m ^i gi ven m 
Table 4 from the condition 

Axo+cfaiJ = Min (mee,E^) [xo+c( m ^, m ee,J2rn Ui )^ - Xo+C,min < 4 . ( 4 -5) 

and equivalently for m ee and Y rn Vi ■ 

The results show that, even for the most restrictive analysis including LSSPS, part 
of the allowed ranges for m„ e in the context of the owCDM + AN rc \ + m u cosmologies 
are within the reach of the KATRIN experiment. On the contrary this is not the case for 
ACDM + mj, models unless only the information of CMB and BAO (or SN) is included. We 
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Figure 4: 95% allowed regions (for 2 dof) in the planes (m Uc ^m u ) and (m ee ,Y^ m f) from the 
global analysis of oscillation data (full regions). We also show superimposed the 95% upper bounds 
on m v from cosmological constraints for the different analysis as labeled in the figure. 



also find that near future neutrinoless double beta decay can test some of the allowed ranges 
in all these scenarios. This will be complementary to the improvement on the expected 
sensitivity from upcoming cosmological probes such as the Planck mission [69]. 

5. Summary 

In this work we have studied the information on the absolute value of the neutrino mass 
which can be obtained from the analysis of the cosmological data in owCDM + AN Te \ + m v 
cosmologies in which besides neutrino masses, one allows for non-vanishing curvature, dark 
energy with equation of state with u ^ — 1 together with the presence of new particle 
physics whose effect on the present cosmological observations can be parametrized in terms 
of additional relativistic degrees of freedom. To break the degeneracies in these models, 
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at least the information from four different cosmological probes must be combined. Thus 
we have performed analysis including the data from CMB experiments, the present day 
Hubble constant HO, measurement, the high-redshift Type-I SN results and the information 
from large scale LSS surveys. We have compared the results from the analysis when the 
full shape information from the LSS matter power spectrum is included versus when only 
the corresponding distance measurement from the baryon acoustic oscillations is consider. 

Our results are summarize in Table 2. Because of the degeneracies present in these 
cosmologies one finds a degradation in the constraints of the standard parameters (ie those 
of the ACDM model) when compared with the analysis performed within the ACDM priors 
for the same set of observables. Concerning the neutrino masses we find that neither of the 
two analysis show any evidence for neutrino mass and the best fit is obtained for J2 m v = 0. 
However the 95 % upper bound obtained when using BAO, ^rriu < 1.5, is tighten by a 
about a factor 2, ^2m u < 0.76, by considering instead the full LSSPS. We have compared 
these results with those obtained from different analysis performed in the framework of 
ACDM+m,, models. The corresponding 95% CL bounds are listed in Table. 3. We find that 
for the same combination of observables CMB+HO+SN+BAO (CMB+HO+SN+LSSPS) 
the bound for a ACDM + m v scenario is ^ m u < 0.61 m v < 0.35) which is a factor ~ 3 
(2) tighter than the corresponding one obtained in owCDM + A7V rel + m v cosmologies. 

Finally we have statistically combined these results with the information on neutrino 
mass differences and mixing from the global analysis of neutrino oscillation experiments 
and we have derived the presently allowed ranges for the two laboratory probes of the 
absolute scale of neutrino mass: the effective neutrino mass in single beta decay m Ue and 
the effective Majorana neutrino mass in neutrinoless /3/3 decay m ee . These results can 
be used to directly address the capabilities of the future j3 and neutrinoless-/?/? decay 
experiments to probe the allowed parameter space. 
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